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Abstract 

We report a measurement of the branching ratio B{D^ (Pt^V()l 
B{D~^ — > K*^(.'^i>() from the Fermilab charm hadroproduction experiment 
E791. Based on signals of 49 it 17 events in the p^e'^Ug mode and 

54 lb 18 events in the — p^p&v^ mode, we measure 

B{D^ p°e+z/e)/Hp+ ^ K*^e^ve) = 0.045 ± 0.014 ± 0.009, and 
B{D+ p^ii+i^^)/B{D+ Z*V+z/^) = 0.051 ± 0.015 ± 0.009. 

Combining the results from both the electronic and muonic modes, we obtain 

B{D^ p°£+Ui)/B{D+ K*^tvi) = 0.047 ± 0.013. 

This result is compared to theoretical predictions. 
13.20.-v, 13.20.Fc, 13.30.Ce 
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Semileptonic charm decays are useful in probing the dynamics of hadronic currents since 
the Cabibbo-Kobayashi-Maskawa matrix elements for the charm sector are well-known from 
unitarity constraints. Form factors for Cabibbo-suppressed (CS) c d semileptonic decays 
can be related via Heavy Quark Effective Theory (HQET) to those for h u semileptonic 
decays at the same four- velocity transfer |T|. Since knowledge of the form factors ixih u 
transitions is vital for extracting Vub from h ^ u semileptonic decays in a model-independent 
way, study oi c ^ d semileptonic decays can improve our knowledge of Vuh- Although con- 
siderable progress has been made in studying CS semileptonic charm decays to pseudoscalar 
mesons 0, the only previous result on CS semileptonic charm decay to a vector meson is 
based on four p^fi^u^ events 0. In this Letter, we report a new measurement from the 

Fermilab hadroproduction experiment E791 of B{D^ p^l^vi^) / B{D^ K*^l'^v^) based 
on more than 100 p^d-^ve. decays in the combined electronic and muonic modes. 

The E791 experiment recorded 2 x 10^° events from 500 GeV/c vr" interactions in 
five thin targets (one platinum, four diamond) separated by gaps of 1.34 to 1.39 cm. Preci- 
sion tracking and vertexing information was provided by 23 silicon microstrip detectors (6 
upstream and 17 downstream of the targets) and 35 drift chamber planes. Momentum was 
measured with two dipole magnets. Two segmented threshold Cerenkov counters provided 
Tt/K separation in the 6 — 60GeV/c momentum range P|. 

Candidates for p^i'^ui, p^ tt+tt" - and D+ K i+ue, K -> K-7i+ decays 

(charge-conjugate states are implied throughout this Letter) are selected by requiring a 
three-prong decay vertex of charge ±1 with one of the decay particles being identified as 
a lepton. A segmented lead and liquid-scintillator calorimeter is used to identify the 
electrons, based on energy deposition and transverse shower shape. The probability that a 
71 (K) is misidentified as an electron is about 0.8% (0.5%). Muon identification is provided 
by two planes of scintillation counters oriented horizontally and vertically, located behind 
shielding with a thickness equivalent to 2.5 meters of steel (15 interaction lengths). All 
the muon candidates are required to have momentum greater than 12 GeV/c to reduce 
background from decays in flight. The probability that a tt {K) is misidentified as a is 
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about 1.6% (2.4%). 

Once the lepton is identified, the other two tracks in the vertex (/ii,/i2) are assigned 
hadron masses. We define the right-sign (RS) sample as vertices in which the lepton 
and candidate have the same charge; hi and /i2 are then oppositely-charged. For 

— > K t^v^ candidates the hadron with odd charge is assigned a kaon mass, while for 
(Pi'^vi candidates h\ and are both assigned pion masses. The RS sample contains 
both signal and backgrounds from reconstruction errors and other charm decay channels. 
The wrong-sign (WS) sample, in which the lepton and candidate have opposite charge 
and h\,h2 have like charge, provides an estimate of the shape of the background under 
the p° peak in the RS sample. Both kaon-assignment hypotheses are kept for WS 

— > K t^vi candidates. 

Due to the undetected neutrino in the decay, there are two solutions for the 
momentum. We choose the lower-momentum solution since Monte Carlo studies show that, 
for the three-prong semileptonic decays, this solution has a slightly larger probability to 
be correct and offers somewhat better momentum resolution. 

To minimize systematic uncertainties, most selection criteria for — > (Pt^v^ are iden- 
tical to those for — > K (.'^vf, exceptions are discussed below. In addition, all criteria 
except those for lepton identification are identical for electronic and muonic decays. The 
common criteria are the following. A decay vertex must be separated from the production 
vertex by at least 20(T;, where ai is the error on the measured separation. The decay vertex 
is required to be at least ba^n outside the nearest solid material, where cr^ is the error on 
the measured distance. The proper decay time for the candidate is required to be less 
than 5 ps. The hadron candidates in the decay are required to have momenta greater than 
6GeV/c. The minimum kinematically-allowed parent mass for the candidate hih2ii'e 
decay, M^in{hih2ii'e) — Pt + \Jpt'^ ^visi is required to lie between 1.6 and 2.0GeV/c^, 
where pr is the transverse momentum of h\h2^ with respect to the flight direction and 
Myis is the invariant mass of hih2i. When masses arc correctly assigned, the M^m distribu- 
tion has a cusp at the mass. This distribution is broadened and shifted to lower mass 
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if there are additional neutral hadrons in the final state. The potential feedthroughs from 
hadronic decays such as D'^ K^tt'^tt'^ and D'^ — > tt'tt'^tt'^ are removed explicitly by ex- 
cluding candidates with either Ktttt or tttttt invariant mass within 20 MeV/c^ of the mass. 
Feedthrough from the Cabibbo-favored (CF) decay Df — > followed by — > K^K~ is 

ehminated by excluding the region between 1.01 and 1.03GeV/c^ in the K'^K~ invariant 
mass. 

In the rarer D'^ — > p^t^vi mode, additional selection criteria are required to reduce non- 
charm background and to eliminate feedthrough from the CF K t^vi mode, which 

has a rate 20 times larger. The following criteria, applied only to — > p^t^v^^ further 
reduce non-charm background. The maximal missing mass squared, M^-^^ — + M^-^ — 
2MD\jMl^g +Pt, is required to be in the range —0.10 to 0.15 (GeV/c^)^. The scalar sum 
of the transverse momenta of the daughter tracks with respect to the flight direction is 
required to be greater than l.OGeV/c. Although these quantities are partially correlated 
with the minimum parent mass, they do provide additional discriminating power. 

When a K from the CF mode — > K*^i~^i>i is misidentified as a tt, the reflected di-pion 
invariant mass is similar in position and shape to the p° resonance. It is thus imperative 
to reduce contamination from — > K i'^ug to a level well below the signal. This is 
achieved with three selection criteria applied to candidate — > p^i^i^i decays, but not 
to the normalizing mode. 1) The minimum parent mass computed for a Kt^Iu^ hypothesis, 
MjniniKniui), is required to be greater than 2.0GeV/c^. Monte Carlo studies show that less 
than 5% of observed K £+z/£ decays populate the Mmin^KTidi'i) distribution above 

2.0GeV/c^, while about 70% of p^t^vi decays populate this region (when a pion is 

incorrectly assigned the mass of a kaon). 2) Information from the Cerenkov counters for both 
hadron candidates is used to reject about 51% of Kt^ pairs, yet keep about 92% of tttt pairs. 
3) Although no significant peak in Kt: invariant mass remains after these requirements, 
the K'K mass for the hadrons is still required to be outside the interval 0.85 to 0.93GeV/c^. 
These three cuts combined with those described earlier result in a relative reduction factor 
of nearly 200 for the K*^ t^vi mode compared with the p^lt^vt mode. 
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Figures |l| and |^ show the signals in p^i^i^e and the normahzing channel 

K i^Ui for both the electronic and muonic modes. Simultaneous binned maxi- 
mum likelihood fits to both the RS and WS distributions are performed separately for the 
electronic and muonic channels in both the p i~^i^£ and K decays. Two 

functions are used in the fit: a p-wave Breit-Wigner shape describes the p° signal and the 
function F{M) = Nq{M — mo)°exp Ci(M — mo) + C2(M — mo)^ , where Nq, mg, a, Ci and 
C2 are free parameters, characterizes the background under the p° (K*^) peak in the RS sam- 
ple which is assumed to have the same distribution as the WS sample. The normalizations 
for the RS background and the WS distribution are allowed to vary independently. In the 
case of —>■ p^i^ue, the shape of the p° is modified by the energy available in the decay. 
Thus, the p° mass and width are taken from Monte Carlo simulation of the p^C-^vi 
decay. For the K i'^U£ mode, both the K width and peak position are free pa- 

rameters in the fit; the values obtained from the fit agree with those from Monte Carlo. 
The p-wave Breit-Wigner functions from the fits are integrated from 0.65 to 0.90 GeV/c^ for 
the vr+TT^ invariant mass for the p^^^i^i signal, and from 0.85 to 0.93GeV/c^ for the 

Kit invariant mass for the K*^i^h'e signal. The yields for both the p^C-'^i'e. and 

^ K l^Vf, channels are listed in Table |. 

The efficiencies are factorized into two parts. The Cerenkov particle identification effi- 
ciencies are determined from a sample of K^-k^-k^ decays from real data, where the 
kaon and pions can be identified by charge alone. The rest of the reconstruction efficiencies 
and acceptances are determined from Monte Carlo simulation. Only the relative efficiencies 
for D+ p^l+Vi and D+ K i+u^ enter our final result. The overall efficiencies for the 
— s> p £'^iy£ and K channels, as well as for the background modes, are listed 
in Table |. 

Only backgrounds which populate the region and mimic the p*^ resonance are trou- 
blesome in the raw p^i~^i^e signal. Using simulated hadronic charm decays from the 
channels most likely to feed into p^i^i^e signal, we found that hadronic charm decay 
feedthrough to the p^i~^i^e signal is negligible. 



The background contributions to the signal mainly come from real semileptonic charm 
decays. The amount of feedthrough from them is based on efficiencies estimated from Monte 
Carlo studies and Particle Data Group (PDG) branching ratios unless otherwise noted. 
To estimate the backgrounds from D'^ decays, the to production cross section ratio 
(Jj^+/(Jd+ is needed. The weighted average of the measurements from hadronic charm pro- 
duction experiments JlOl with a conservative error, 0.58 ±0.15, is used. The most significant 



semileptonic charm decay backgrounds are listed in Table H along with the correspond- 
ing estimated number of events in the p^l^vi signal. Requiring candidates to lie 
between 0.65 and 0.90GeV/c^ in vr+vr" mass effectively removes decays from — > (pi^ue, 
(f) K^K~ , as well as from V^^^e^ ~^ V^^'^t- The contribution to the signal from 
these modes is negligible. The backgrounds from non-resonant or higher-mass-resonance de- 
cays are negligible as well. 

After background subtraction, the final numbers of signal events are 49 ± 17 for 
—>■ p^e^Vf. and 54 ± 18 for p^fi^u^. The yields in the normalizing channels are 

892 ± 52 for D+ Z*°e+z/e and 769 ± 54 for D+ K*°fi+u^. 

Systematic errors associated with lepton identification are largely cancelled in the ra- 
tio of the p i~^i^e and — K i'^i'£ decay rates. Remaining sources of systematic 

error are 1) uncertainties in the branching ratios used in background subtraction, 2) un- 
certainty in the to production cross section ratio a^+/a£i+, 3) determination of 
relative efficiencies and 4) the fitting procedure. The uncertainties in the relative efficiencies 
are dominated by the momentum dependence of the Cerenkov identification and the depen- 
dence of the — s> p^t^vi acceptances on the form factors. The effects of the uncertainties 
in the assumptions made in the fit are evaluated by varying the width of the p° peak and 
the shape of the background distribution. For the electronic channel, the four sources con- 
tribute approximately equally to the systematic error, each with a size about one third of the 
statistical error. For the muonic channel, the first three sources contribute approximately 
equally (about one third of the statistical error each) while the last source contributes an 
uncertainty about one sixth that of the statistical error. Further studies were performed to 



search for other potential systematic effects by varying the selection criteria one at a time. 
No significant effects were found. 

The rate for the decay D Vi^Ui, where \^ is a vector meson, is determined in the 
limit of massless i by three form factors Ai{q'^), A2{q'^) and V{q'^), where is the square 
of the four-momentum transfer from D to V The present experimental information 



for — s> K i'^Ui decays is usually presented as Ai(0), and the ratios of form factors 
r2(0) = A2{0) / Ai{0) and ry(0) = V (0) / Ai{0) , with an assumed dependence proportional 
to (1 — q'^/Mp)"^ and Mp ~ 2.1 to 2.5 GeV/c^. Since we have insufficient statistics in 
this experiment to measure the — > p^l^vi form factors, the simulations on which our 
efficiencies are based assumed r2(0) = 0.82 and ry(0) = 2.0, close to the present world 
averages for K t^vi. We have checked the effect on the detection efficiency of 

significantly different values of the input and vy. Specifically, assuming r2(0) = 0.0 or 
ry(0) = 1.0 changes our p^t^vi branching fraction by less than (10 ± 5)% of itself. 

These variations have been included in the systematic uncertainty. 

From the background-subtracted event yields and the efficiencies for (Pt^vs^ and 

^ K i^ui decays, the following branching ratios are determined: 

(4.5 ± 1.4 ±0.9)%, 
(5.1 ± 1.5 ±0.9)%. 

We combine the results from the electronic and muonic modes, taking correlated errors into 
account, to obtain a final result of 

where the error includes both statistical and systematic uncertainties. In Table ^ our result 
is compared with the only previously-published experimental result, from E653 and 
various theoretical predictions. Our result, which is sensitive mainly to the form factor 
Ai at q^ ~ 0.5 (GeV/c)^, agrees only marginally with quark model predictions [^,|T^,|T3 
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but agrees well with recent lattice QCD calculations [114Hl7|] and several other theoretical 
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predictions [|T8H20[]. It thus begins to discriminate among models that are also used to 
predict form factors for 6 ^ m semileptonic decays to extract V^b- 
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TABLES 

TABLE L Numbers of fPt^vn and — > K"'^ i^V(^ signal events from the fit to the 

data, estimated numbers of background events, efficiencies E for each decay mode and input to the 
calculation of backgrounds. Efficiencies include branching ratios for non-charm decays. 
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Kaw p 


lUi Signal 
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0.005 


ISGW2 


^ prediction for 
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0.004 
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D+ ^ 
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54 lb 18 


0.19 


Background subtracted signal 
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49 lb 17 


0.16 







Normalizing mode 



D+ K e+vc n 769 lb 54 0.19 Fit results for 

K~7T~^ e 892 lb 52 0.19 normalizing signals 
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TABLE II. Comparison of our results with the only previously-published experimental re- 
sult and theoretical predictions. The PDG branching ratio for — > K*^£'^vi and 
lifetime are used to calculate the experimental decay rate for p^l'^v^. Most theoretical 

results are calculated for decays. To compare these results for decay with the experi- 
mental results for D'^ decay we use the relations T{D+ K i+ui) = r(Z)° ^ K*-e+Ui) and 
T[D~^ — > p^t^Vf) = 1/2 X T{D^ p^t^Vf), where the factor of 2 difference between the two decay 
rates arises from the l/\/2 coupling of to the p^ . The second column indicates the method used 
to obtain the results, where QM stands for quark model, HQET for heavy quark effective theory, 
SR for QCD sum rule, and LQCD for lattice QCD. 



Group 
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I 




T{D+ ^ p^£+ue) (IQiOs-i) 


E791(this work) 


Exp. 


e 


0.045 ± 0.014 ± 0.009 


0.20 ±0.07 ±0.05 


E791(this work) 


Exp. 




0.051 ±0.015 ±0.009 


0.22 ±0.07 ±0.05 


E791(this work) 


Exp. 




0.047 ±0.013 


0.21 ±0.06 
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QM 




0.022 


0.12 
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0.06 


0.225 


Ball 1^ 


SR 
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FIG. 1. Mtttt distribution for D+ p'^i'^i^e candidates. The vertical arrows indicate the mass 
window for the final p^£^ui candidates. For each leptonic mode, a simultaneous fit is made 

to the RS and WS data. The shape for the background distribution in the RS data is constrained 
to be the same as that of WS distribution, but the relative normalization is allowed to vary. 
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FIG. 2. Mktv distribution for K i'^u^ candidates. The vertical arrows indicate the 

mass window for the final I?+ — ^ if candidates. For each leptonic mode, a simultaneous fit 

is made to the RS and WS data. The shape for the background distribution in the RS data is 
constrained to be the same as that of WS distribution, but the relative normalization is allowed to 
vary. 
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